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ABSTRACT 

We solve general 1-matrix models without taking the double scaling 
limit. A method of computing generating functions is presented. We 
calculate the generating functions for a simple and double torus. Our 
method is also applicable to more higher genus. Each generating function 
can be expressed by a "specific heat" function for sphere. Universal 
terms, which survived in the double scaling limit can be easily picked 
out from our exact solutions. We also find that the regular part of the 
spherical generating function is at most bilinear in coupling constants of 
source terms. 
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1 Introduction 



Matrix models JTJ 0, Q ||] in the double scaling limit |IJ enable us to treat two 
dimensional gravity nonperturbatively. This special limit is necessary for us to 
investigate the universality of matrix models. Many aspects of universality hidden 
in these simple models have been revealed by a lot of people. Matrix models in the 
double scaling limit are described by the string equations, a nonlinear differential 
equation for the full specific heat function. The topological expansion of this specific 
heat requires a large cosmological constant t. If one want to define the specific heat 
for small t, nonperturbative input must be demanded. This point makes the physical 
meaning of the double scaling limit unclear. 

Our aim of this article is to solve general 1-matrix models exactly as random 
lattice models. Until now matrix models have been solved in the double scaling 
limit or in simple potential cases. Our solution of general 1-matrix models are exact 
even when the matrix size N is finite. All of a 1-matrix model can be shut in a 
function. We find the generating function for each topology when the function is 
given. To study non-universal aspects of 1-matrix models, we rearrange the notion 
of the scaling operators || in these discrete models. Using this refined formalism, 
we can understand a lot about the bulk contributions and the finite size effects 
of general 1-matrix models. Our exact solution for general 1-matrix models are 
important, since these are precious examples of exactly solvable lattice models as 
well as toy models of gravity. We can easily pick out contributions that survive in the 
double scaling limit from our solutions. This approach to the universal contributions 
gives us an explicit answer to the leading generating functions. As we can include 
source terms of the physical observables, we can understand the correlation function 
of these observables more directly than the double scaling limit approach. We cannot 
find equations like the string equations in our formalism, however. We may rather 
consider that the double scaling limit approach and ours are complements each of 
the other. 

Before we show main results, we briefly review known results given by Bessis, 
Itzykson and Zuber (BIZ) 0. A general 1-matrix model is defined by a integral 
over a N x N hermitian matrix $, 



We always fix the coefficient of tr$ 2 to 1/2 in the action. If we rearrange and fine 
tune the coupling coefficients {A}, the system behaves as one of the Kazakov series 
with some source terms. 
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Using the orthogonal polynomial technique, the generating function of this model 

is0 

j? rs\\\ 1 i Z "({A}) 
En({\}) = ^log 



N 2 to Z N (Q) 

-yd--) iog W) + iiogMffl (i 2) 



If we denote the orthogonal polynomial of degree k as Pk(fJt), i?fc({A}) is defined by 
a recursion relation for the polynomials, 

fiP k (fi) = P k+ M + R h P k -M ■ (1.3) 

/i in ( |1.2| ) is defined as, 

/oo 
(1.4) 

where 

The most characteristic point of this model is that the large- N expansion of is 
a genus expansion 

oo 

E N = N~ 2h e {h) , (1.6) 

h=0 

where e^ h > is a contribution from genus h. 

The authors of || expanded £jv using the Euler-Maclaurin formula and obtained 
e C0 formally. Let us introduce some notations to show their result compactly. It is 
convenient to note e = 1/N, x = k/N and introduce r e and its e-expansion, 

oo 

r e (x, {A}) = R k ({X})/N = £ e 2s r 2s (x, {A}) . (1.7) 
r e satisfies a recursion relation derived from (|L1 

oo 

x = r e (x,{A}) l + ]T2(p + l)A p J2 r e {x + Sl e,{X})---r e (x + s p e,{X}) 
p =1 paths 

' (1-8) 

where Spaths means a sum over the 2 P +iC P paths along a staircase depicted in Fig.l 
and integers Sj, i = 1,2, ...,p are the height of the stair when i-th descending 
down. For example, we can associate {s\, s%, S3, S4, S5} = { — 1, —1, —2, 0, 1} for the 
path drawn in Fig.l. 



1 Our definition for En is slightly different from Eq.(4.3) of M. Their sign of e^ 1 ' is also opposite 
to ours. Take care when comparing their results with ours. 
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Fig.l: An example of path from point A to point B. There are 2p+iC p paths 
which connect two diagonal points of a lattice of size p X (p + 1). In this 
example, p = 5. We assign a set of p integers for each of the paths. 

Expanding both sides of this recursion relation in e 2 and comparing the coefficients 
of e 2s , we can solve r 2s . If we take a limit e — > in ( |1.8|) , we get a function W(r ) 
which have all information of the model, 



x = r Q 



1 , ^ (2p + 2 r , „ 



V' 



Note that 



W(0) = 0, W'(0) = 1, lim 



= W(r ) 
W(r ) _ 1 ^ 



1.9) 



One of our goal is to understand for given W{tq). We introduce two new 
notations. The first is defined by a large- N expansion 



log 



r e (x, {A}) 
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x 



E^2l f 2 S (a;,{A}) 



We show the first few of them: 
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In our formalism, these quantities play an important role. The next is associated 
with the last term of the right hand side of (|1.2| ) p|, 



1 

2N 



' h ({X},N) h ({X},-N) 
lo S i 77T\ log ■ 



ho(0) 



ho(0) 



s=l 



;i.i2) 



Since A p is rewritten by a derivative of W(r ) at r = 0, 



h ({\},N) 



dp, e 



-V(j*,N) 



-oo 
oo 



dfx exp 



P < ty»+"(o) 

h 1 



{(2p + 2)\ Np 



(1.13) 



The first few of them are: 

K = 0, K 2 = ~W"(0), 

K A = -— \22{W"(0)} 3 - 20W"(0)W"\0) + 3W {4) (0) 
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Then the contribution from genus h is written as 



;i.i4) 



o 



dx(l-x)f 2h (x,{\}) + K 2h ({\}) 



h 



-E(-ir 



(2«)! 

where are the Bernoulli numbers 



{(l-x)r 2 (hr-,)(x,{\})} 



(2.-1) 



;i.l5) 



-Bi = — , i?2 = — i 5'i = — ..... 

6 30' 42' 



;i.l6) 



In ( 1.15|) the superscript (2s — 1) stands for the (2s — l)-th derivative with respect to 
x and the symbol |J means taking the difference between the values of the function 
at x = 1 and x = 0. We can directly solve the recursion relation for r 2s and 
calculate the generating functions from (|1.15 ) when W{tq) is simple. It is, however, 
a troublesome task to find directly from (|1.15| ) even when W(r ) is quadratic. 
There are roughly two way to advance: to take the double scaling limit and drop all 
non-universal information, or to solve ( |1 . 1 5| ) for general cases more explicitly. Our 
method achieves the latter. 

The main outcome of this article is the generating functions for a simple and 
double torus. The answer for sphere has already been given in 0]. If we use a 2 such 
that W(a 2 ) = 1, the answer for is very simple, 



1 

12 



\og{a 2 W'(a 2 )}. 



;i.i7) 
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On the other hand is a little complex, 



,(2) 



a 2 



28{T^"(a 2 )} 3 - 29W'(a 2 )W"(a 2 )W"'(a 2 ) + 5{W'(a 2 )} 2 W^(a 2 



UA0{W'(a 2 )} 5 
9{W"(a 2 )} 2 - 4W'(a 2 )W'"(a 2 ) 
2880{W'(a 2 )} 4 
1 1 W"(a 2 ) 1 1 



240 240 a 2 {W'{a 2 )Y 240 a 4 {W'(a 2 )} 2 



Though we do not compute e^ h ' for h > 3 this time, there are no obstacles. They 
will be more complex than . These explicit answers for e^' allows us to compute 
exact correlation functions of scaling operators. 

The outline of this paper is as follows. In section 2 we make a systematic study 
on the representations of scaling operators by the matrix $. We write down the 
polynomial W(r) in which all information about the k-th model with source terms 
are packed. The formulation in this section shows its ability when we study non- 
universal aspects of matrix models as well as the universality. In section 3 we 
compute the spherical generating function for the general setting of section 2 and 
show that the regular part is naturally separated from the result. Section 4 is the 
main part of this paper. We give here a method to perform the program of BIZ 
thoroughly for general W(r). It is essential to understand weighted summations 
over paths appeared in the formulation of BIZ. We solve them using generating 
functions for these sums. As an application of the result of section 4, we compute 
the correlation functions of scaling operators for a simple and double torus in section 
5. We also study the leading singular terms which survive in the double scaling limit. 
We discuss advantages and disadvantages of our method in section 6. In appendix 
A we prove that a generating function of a special weighted sums over paths is a 
g-deformed combination. Appendix B is devoted to another generating function for 
sums over the paths. 



2 Fundamental Representations for Scaling Op- 
erators 

The scaling operators in 1-matrix models have simple scaling properties This 
scaling properties are specified by the string susceptibility 7 and gravitational di- 
mensions of scaling operators. The k-th. criticality is realized by a fine tuning of the 
parameters {A} in The string susceptibility for the k-th model is —1/k. It 

is well-known that the gravitational dimension of the scaling operator in the k-th 
model is a multiple of the absolute value of the string susceptibility I7I —1/k. We 
can express these objects by the matrix $. There are some ambiguities in how to 
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represent these scaling operators with respect to $. Generally, a scaling operator 
with gravitational dimension d is written in a form of a linear combination of the 
following infinite set of "fundamental representations" , 

p%[% d=^ n,M = 0,1,2,... . (2.1) 

The definition of these fundamental representations is simply given by a generation 
rule 

= /$-!)/*[*] - PJSW*] . (2-2) 

„«[*] = (jl±iHj£±^^ t r^. (2.3) 
Ho 1 J (2M + 4)! v ' 



This generation rule (2.2) means that each fundamental representation is a linear 
combination of puncture operators (|2.3|) 



PnM$] = E( „ )(-i) p Po +p m, (2-4) 

where 



n \ _ n! 



\ P J p\(n — p)\ 

The number of puncture operators in (|2.4j) corresponds to the gravitational dimen- 
sion d = n/k of the scaling operator. The index M specifies the puncture operator 
with the lowest degree when expanded by the puncture operators. In Section 5 we 
will show that the scaling operators with different M degenerate in the double scal- 
ing limit. This means that the index M is non-universal. For example, p\ in pure 
gravity (k = 2) is 

p\m = P im-2 P im+pim 

1 ■ tr $ 6 — tr $ 8 H 1 — tr $ 10 



480 2240 40320 

= — — tr f $ 10 - 18$ 8 + 84$ 6N ) . (2.5) 
40320 v J K J 

Using these fundamental representations, we can formulate a general model, 
which is fine tuned in the fc-th critical point and has (m + 1) source terms. We 
reserve the 0-th entry as the cosmological constant term. The action for this model 
is 

m 

Sm=S a [*]+ y £T if % /h [*] t (2.6) 

where the fine tuned part Scr[^] i s formally written by a special scaling operator 
with negative M and dimension one, 

S cr [$] = iV- pr 2 [$] 

tT V _ y> fk \ l 1)p H!p1^ (2 7) 
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The existence of this special scaling operator prevents us define fundamental repre- 
sentations with negative M. These irregular representations change the coefficient 
of tr$ 2 in ( |2.6| ). We can avoid these representations by redefinition of $ and the 
coupling parameters Tj. If we formally define that n_i = k, M_i = —2, r_i = — 1, 
( |2.6| ) becomes simpler, 

m 

Sm=N+ £ r iP ^M- (2-8) 
t=-l 

We can translate this action in the form of W(r) appeared in the previous section, 

m / r \ 2+Mi / r \ri i 

^)=i-Ey (2.9) 

This fact justifies our definition of the fundamental representations. The definition 
( |2.2j ) has a simple meaning: each scaling operator is constructed by cancelling the 
leading singularity of the two distinct scaling operators with dimension subtracted 
by |7| = 1/k from the original one. Our setting of the 1-matrix models given by 
(|2.9| ) is very useful in later analysis of generating functions. 

Strictly speaking, the definition for the fundamental representations depends on 
the choice of the fine tuned part of the action 5* cr . There is a degree of freedom 
for deformation of the fundamental representation p^/k- This degree of freedom is 
related to the choice of S cr . The modified generation rule is similar to the previous 
one 

/>"*[*, &] = P? n - m l*,b] - pgtl^.b] , 

,o"[<M = ( ^fff 4) t 2) ' (* + »- M -'tr^. (2,0) 

The deformation parameter b must be greater than —k. When 6 = 0, fl2.10|) gets 
back to the generation rule ( j2.2|) . These deformed representations correspond to a 
fine tuned action 

S cr [$, b]=N- pr 2 [$, b] + bp^[<5> , b] . (2.11) 

This additive term is a special marginal scaling operators as before. We assign 
M_ 2 = —1, n_2 = k, r_ 2 = b, for the additive term in (|2.11| ). Then the whole 
action is translated to 



1 k + b)\ k + bj f^ \k + b) \ k + b 

r \ Mi+2 



T; 



-2 



We can easily see that p^/ki^i 1] = Pn/(k+i) 0]- Therefore the parameter b connects 
the k-th criticality with the (k + l)-th criticality of 1-matrix models. We can add 
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ordinary marginal operators, i.e. , operators with gravitational dimension one and 
non-negative M to the critical action S cv . These terms, however, do not cause 
further modification in the fundamental representations. Since a coupling constant 
Tj do not scale when its dimension is one, we should regard these marginal operators 
as source terms. In this point of view, the most general fine tuned action consists of 
marginal operators with M = —2 and M = — 1. The special marginal operator p± 2 
is always necessary to cancel N in the fine tuned action. In section 5, we will show 
that this parameter b (if b ^ 1) is non- universal as the index M. From now on, we 
will mainly discuss the case of b = for simplicity. 

We conclude this section by classifying quantities which characterize general 1- 
matrix models. 

1. Universal quantities. 

k specifies the criticality of the system. The string susceptibility 7 

depends only on this integer, 7 = — 1/k. 

rii is related to a gravitational dimension of a scaling operator, d$ = 

Tli/k. 

2. Non-universal quantities. 

b modifies the fine tuned part of the action, b also connects adjacent 

criticalities. 

Mi shows degrees of freedom in representing scaling operators with the 

matrix $. A scaling operator is constructed by taking a difference 
of the two scaling operators with adjacent values of Mj and with 
gravitational dimension lower than the original one by 1/k. 



3 Non-universal Part in Generating Function for 
Sphere 

The planar contribution of the generating function has been argued since early times 
H], D, ||]. In this section we represent the scaling nature of spherical part without 
introducing any scaling units such as the cosmological constant. 

The spherical generation function is written by the function W(r) as 

e (°) = [ 1 dx(l-x)\og- 
Jo X 

= - drW\r)[l-W{r)}\og\^^\ . (3.1) 
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Integrating by parts, we can rewrite this in the following form, 



e (°) = ^loga 2 + ^ 



2 Jo 



where W is 



W(r) = 



drrW'(r)[2-W(r)} 
W(r) 



If we substitute the most general form for W(r 

W(r] 



, >fi (2p + 2)\ p+1 



(3.2) 



(3.3) 



(3.4) 



to ( p.2|) and eliminate Ai using W(a 2 ) = 1, we have 



,(o) 



iloga 2 + l(l-a 2 )(9-a 2 ) 



+ E 

p=2 



(p-l)(p + 5)(2p + l)! 



a2p+2 _ (p-l)(p + 6)(2p + 2)! ^ +4 



A, 



2(p + 3)(p + l)!p! 12(p + 3)!p! 

, v V (p-l)(g-l)(p + g + 6)(2p + 2)!(2g + 2)! 2p+2q+ 
hhs(p + 3)(q + 3)(p + q + 2)(p+l)lpl(q+l)\q\ p q -\ ■ ) 



p=2 <j=2 



Here Ai is nothing but the cosmological constant. When A p = 0, (p = 2, 3, . . .), this 
result is the same as that of 0. This form, however, is not fit for understanding 
the scaling property of the generating function. We must fine tune to the k-th 
criticality and introduce source terms to understand universality hidden in matrix 
models. The setting of the previous section, 



W{r) 



1 



z + E X 



i=0 



Mi+2 



1 ~ I 



(3.6) 



is suitable for our aim. Substituting this function in ( |3.2| ), we visualize the scaling 
behavior of e^ -* roughly. The integrand in Q3.2|) is described by powers of 1 — r / k. 
The {r}-independent terms in e^ ) are 



1, 2 1 

-log a +- 
2 2 Jo 



dr 



k-l 



k 



E 



p=0 



1 J* 1 



- log fc H ^) 

2 S 4 2±-<p 



a 



2\ k 



p=l 
1 



4 



1 



a 
1 



k) 



-i oo 

-E — 

2 p=0 2 ^ 
2\ 2fc 



fc-1 



P 



a 
1 



r 

~ k 

2 \ 2fc+l+p 



(3.7) 



Picking up linear terms in {r}, we find 



rn 
i=0 



1 



1 - 



fc-i' 



Mi+2 
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rh ( _ r 
k V k 



m 
i=0 



y - 

^ k + m + l + p 



+ 



1 / a 



2 \ M,+2 



2 V k 



1 



a 



k 



+ 



P 

1 fa 2 



i 2 \ k+rii+l+p 

1 ~ t) 



- 1 



a 



2 \ fc+rij' 



2 V fc 



To derive the second and the third term of the most left hand side of 
a property of binomial coefficients such as 



(3.f 



we used 



Mi + 1 



(Mi + m + 2) 



p I v " ' v p — 1 

The bilinear terms in {r} are 

TiTj j n dr( 

i=0 j=0 

X 



Hi 



Mi + 2 
P 



(Mi + n. + 2-p). (3.9) 



r \ " J Mj ■ ii j ■ 1 / r \ " // 
1 ~ k) I fc V 1 ~ fc 



EE r ^ 



i=0 j =0 



1 



Afj+Mj+3 



•1? 



i 1-- 

fc V k 



Mi + Mj + 3 



1 ~ t) 



q n, + rij + 1 + p 
- 1 



P 



1 /a 2 \ Ml+Mj+4 / a 2 \ ni+n ' 



(3.10) 



Here we again used the above nature of binomial coefficients. Introducing the fol- 
lowing rules for binomial coefficients 



1 



P 



= (-l) p (forp>0), ( p J=° (forp>n), (3.11) 
and using W(a 2 ) = 1, we can assemble all contributions in a compact form, 



1. . 3 1 ~ 1 

- log k H > - 

2 4 2 p ttp 



1 J 



m, m 



-o E E E 



-i? 



X 



-ii=-i 

2 \ n.i+nj + 1+p 



/ Mi + Mj + 3 
p= o n» + rij + 1 + p \ P 



1 - 



a 



- 1 



(3.12) 



where negatively labeled quantities This final result naturally separates into the 
regular part and the rest, 

AO) _ JO) , JO) /oioN 

e ~ e bulk + e sing ' ^- L6 > 
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where 




2t ^ o ij ^ n t + n J + l+p{ V 
1 3 ™ (Mi + l)\(ni + k)\ 



Mi + 1 



) 




(3.14) 



and 




o E E T i r i E 



r= — 1 jf= — 1 p=0 



mm oo 



Hi + Uj + 1 + p 



( 



Mi + M 3 - + 3 
P 



) 




2 \ ni+rij+l+p 



(3.15) 



The regular part e^^jjj has a remarkable nature. It is at most bilinear in the coupling 
constants {r} and all the coupling constants are equal in ej^^.. This result for the 
non-universal terms is nontrivial, for we apt to think that the structure of the non- 
universal part changes its form according to the operator content of the model. This 
naive estimate is not true, however, and there is a universal structure even in the 
non-universal part. We should consider this universal structure in the bulk part 
as one of the characteristic of the matrix models. Though we have no proof, it is 
natural to think that there is the same structure in the non-universal part even in 
2-matrix models. On the other hand we regard as the singular part, because 
the factor 1 — a 2 /k becomes singular with respective to a scaling unit which will be 
introduced later. Though we have not yet introduced such a scaling unit, we can see 
a sign of scaling nature through the factor 1 — a 2 /k. We should note that there is a 
simple relation between ej^^ and e^g. If we replace a 2 with in , it equals 

to — e bulk U P ^° { r Independent constant (infinite) value. 

Our result for the bulk part is crucial in matrix models modified by trace-squared 
terms @, 0, || . The effective action for baby universes associated with scaling oper- 
ators depends deeply on this bulk structure of unmodified models. Our consequence 
shows that the Liouville dressing {3 for a source term always changes its branch 
when the coupling constant of a trace-squared term associated with the source is 
fine tuned. This also denies the existence of an upper bound for the gravitational 
dimension of scaling operators |Q. 

When we deform the system by b as in the section 2, the generating function is 
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compactly given in the following form 



1 3 1 00 1 

e(°)(6) = i]og(* + 6) + 



p=i 



m m 



-o E E 



X 



i=-2j=-2 p 

2 \ n.i+rij+1+p 

? / 

1 



a 



k + b 



J Q ni + nj + l+p 
1 



Mi + Mj + 3 
p 



(3.16) 



where we used notations of section 2 for negatively labeled quantities. It is important 
that the structure of the spherical generating function (|3.13|) is kept even when 
nonzero b. One can easily check that when b = 1, the system behaves as the (k + 1)- 
th Kazakov series, 



c W(6 = i) = ±iog(* + i) + £-±£;± 



i ^ i 



X ^ 1) P 



2 TiTi S rii + rij + 1 + p 



i=-ij=-i 



1 - 



p=0 
rii+rij+l+p 



Mi + Mj + 3 
P 



fc + 1 



- 1 



(3.17) 



where n_i = k + 1 instead of k. 



4 Generating Functions for Higher Genus 

In this section we compute generating functions for higher genus with respective to 
the function W(r) in which all the information of the model are packed. We must 
overcome a difficulty in expressing the quantity r 2s appeared in in terms of W(r ) 
to accomplish this task. We realize this with the aid of a g-deformed combination 
and a generating function of summations over paths connecting two fixed points 
(See Fig.l in §1). 



4.1 Torus 



Our starting point for the toric contribution of the generating function is read from 



,(1) 



r dx{1 - x) r .i- V'(o) +i {a - X ) log ^ j 

Jo r 8 12 I x J 



(i) 



(4.1) 



where x is nothing but the function W(tq). The last term of fl4.1|) can be written 
only by a 2 and the derivatives of W(r Q ) at r = 0. To show this we define derivatives 
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of r (x) with respect to X as rj^ = (d n /dx n )r . If we differentiate both sides of 
(|1.9Q with respect to x, we have 



rf* = [W'iro)}- 1 . 



(4.2) 



So we have 



1 -a;) log 



r 



x 



(i) 



r 1 - x (i) 1 - x 
- log 1 rk 

x r x 

, 1 - W(r ) 1 - W(r ) 

g W(r ) ro^'(ro) W(r„) ' 1 ' ' 



Since jy(a 2 ) = 1, the difference between the values of the above function at x — 1 
and x = is 



1 

12 



1 - x) log 



(1) 



X 



— log a 2 lim 

12 & 12 r ^0 



1 - W(r ) 1 - W(r ) 



r W'(r ) W(r ) 
Using the facts like W(0) = 0, W(0) = 1, we can find the limit value, 



lim 

ro— >0 



1 - W{r ) 1 - W(r ) 
r W'(r ) W(r ) 



lim 

ro— >0 



lim 



1 



1 



r W'(r ) W(r ) - W(0) 



ro-o {r W'(0)} 2 V 1 + 2 
W"(0) 



(4.4) 



(4.5) 



Thus we have 



— ((1 — x) log — 1 
12 l v ; xj 



(i) 



-— loga 2 + — W"(0). 
12 & 24 v ; 



(4.6) 



We now rewrite r2/ro in terms of W(r$). To do this we solve the recursion 
relation for r e (|1.8Q . We identify the coefficient of e 2 on both sides of ( |i~8|) , and find 







r 2 + ]T2(p+l)A p 



2p+ 1 
P 



r 2 ro + r 



+ E E^ 



s 2 rj 2) rg x + 



1 



paths 



7i E E' 

paths l J 



2p+ 1 
P 

■?>)V 2 



p-i 



(4.7) 



where Y^ij ' means that the case i — j is excluded in the sum. There are two unknown 
summations over paths. We must find their values for given p. It is difficult to guess 
the values of these sums for general p from the individual values of the sums at small 
p. We need some clever devices to understand these sums. 
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Let us introduce a general definition for such sums over the paths. For a given 
partition of an even integer 2n, 



2n = 1 • at + 2 • a 2 H h (2n) • a 2n , 

we define a summation over the paths as 





' 2n 


-1 


" 2n 


I p (l ai ■ 2 a2 ■ ■ ■ (2n) a2n ) = 




E 






.1 =1 


paths 


.i =i 



El n\ 



, n 2 
k 2 



(4.9) 



In (f4.9|), ' means a summation over all possible m integers 1 < ki, k 2 , ■ ■ ■ , k m < p 
such that any two of them do not coincide, m integers 1 < n\ < n 2 < - ■ ■ < n m < 2n 
are defined so that ai of them are 1, a 2 of them are 2 and so on. Then we have some 
relations 

2n 2n m 

^a g = m, ^qag = ^2fii = 2n . (4.10) 

q=l (?=1 i=l 

These summation vanish in the case when 2n is a odd integer, because I p should 
not change its value when we change the sign of the height — * — sj^. The two 
summations appeared in ( |4.7|) are simply described as 



Jp(2) 



1 p 
paths 4=1 



.) E E s * s j ' 

paths *ii 



(4.11) 



where we omit entries whose a« are 0. Thus all we have to do is describe I p {2) and 
ip(l 2 ) for general p. 

Two generating functions give us a key to solve this problem. To our surprise, one 
of these generating functions is a g-deformed combination or a g-deformed binomial 
coefficient. The definition of it is similar to an ordinary combination 



p 2p+l 

1] (1 + e"y) = £ 2 P+ iC r (xW 

s=— p r=0 



(4.12) 



where p is a non- negative half integer. It is easy to check that this function 2p+ iC r (x) 
is obtained by replacing integers with g-deformed integers (q = exp(x/2)) in the 
usual expression for binomial coefficients, 



2p+l 



C r (x) 



N 2p+1 (x)\ 



N r (x)\N : 



2p+l-r 



»! 



(4.13) 



,,, , . , rx . x 
JM r (x) = sinn — / smn — , 
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N r {x)\ = Nx(x) ■ N 2 {x) ■ ■ ■ N r {x) . 
We give some characteristics of the g-deformed combination, 



2p+lC*r(0) — 2p+\C r , 
2p+l p 



Y J 2 P+ xC r {x) = 2^ I] cosh^, 

r=0 r=—p 

2p+l 

J2(-iy 2p+1 c r (x) = o. 

r=0 

As explained in appendix A, 2 P +iC p (x) is nothing but a generating function of a 
following summation over the paths, 



e n 

paths 4=1 



e StX = 2p+iCJx) 



N 2 p+i(x)\ 
N p (x)\N p+1 (x)\ ' 



(4.14) 



Now we call the expansion coefficient of x in the g-deformed combination G p (2n) 



2p 



p+\C p {x) — ^2 G p (2 



n)x 



2n 



n=0 



This coefficient G p (2n) is the sum of I p over all partitions of 2n, 

1 i p ^ 2n 

G p (2n) ~ ~ 



paths 
all partitions of 2n 



(4.15) 



1 / p 

E 2 lE s * 

paths " =1 



Thus G p (2) is the sum of I p {2) and J p (l 2 

G p (2) = 

paths 
= J p (2) + / P (l 2 ). 

Since Taylor expansion of a g-deformed integer is 

AUx) = r fl + — (r 2 - l)x 2 + — — (r 2 - l)(3r 2 - 7)x 4 + 0(x e 
r\ ) [_ -r 24 v ; 576Q v a ; v 



(4.16) 



(4.17) 



we can compute the summation over paths in ( 4.16 ) 

/2p+l p p+lN 



G P {2) 



2p + 1 W 

P / \ i , 



IE-E-E 



12 



(P + 1) 2 P 



2p + l 
P 



(4.18) 
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We can check this result for small p. 

As we find the sum I p {2) + / p (l 2 ), we must find one of them, say I p {2). To do 
this we define a generating function of I p (2n) 

r„(*) = E E e Si!E 



paths 



i=l 



= ]T/ p (2n)x 2n . (4.19) 

n=0 

There are no odd terms in x due to the symmetry Sj — Sj. In appendix B we show 
that is given by a summation of g-deformed integers weighted by a binomial 

coefficient 

' "~ ' / p+l+i-j ) ( p-l-i + j \ e ( P -i-j) x 



r P (x) = EE 

i=0 j=0 

p-1 



3 



E ( \ + 1 J N 2p+l . 2n (x) . (4.20) 



n=0 



The proof of this relation is interesting. It is hard to understand this relation alge- 
braically. We prove it in appendix B diagrammatically. Comparing the coefficient 
of x 2 in both sides of ( |4.20| ), we find I p {2) (see appendix B for details of calculation), 

I p (2) = ^EV 2 ^ 1 y2p+l-2n){(2p+l-2n) 2 -l} 

= ^(P+ 1 )p( 2P p 1 ) • (4 ' 21) 
Thus the rest I P (1 2 ) is 

J p (l 2 ) = G p (2) - I p (2) = ^(P+ l)p(p ~ 1) ( 1 ) • (4.22) 

This vanishes when p = 1 as expected. 

Since both of I p {2) and / P (l 2 ) have a factor (p + l)p and a combination 2 P +iC p , 
we can rewrite the recursion relation for r 2 ( |4.7| ) with respect to derivatives of W(r ) 

r 2 W'(r ) + l4 2) r W"(r ) + ^{r^ fr^W'"^) = . (4.23) 
o 1/ 

We already know that = [W( r o)] _1 - As for 7q , we rewrite it with derivatives 
of W(r ) differentiating both sides of x — W(r ) twice with respective to x, 

rj 2) = -W"(r )[W'(r )]- 3 . (4.24) 
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At last we succeed in describing r 2 /r by W{ro) 



12 
r 



W'(r ) 
1 / 



l{W"(r )} 2 1 W'"(r ) 



6{W'(r )} 3 12{W'(r )} 2 

2 

logW(ro)- 



d 



12 \W'(r )dr J vu/ - (4 ' 25) 

Thus we can explicitly perform the integration in the first term of the right hand 
side of (|OD 

fdx{l-xp = -^r\r{l-W{T))^-^^-\ogW'{r) 
Jo r 12 Jo drW'(r)dr 



= —W"(0)-—\ogW'(a 2 ) 
12 v ; 12 & v ' 

Collecting all contributions, we get to the final answer 



(4.26) 



^"(0) -i log W'(a 2 )-±W"(0) 



1 

12 



\og{a 2 W'(a 2 )}. 



(4.27) 



The result is surprisingly simple. A matter of special mention here is cancellation 
of terms which include W" (0). We can easily check our result by a simple example 
of § 

W{r) = r + 12Air 2 , (4.28) 
where Ai is a coefficient of tr $ 4 in the action. Since 



a 2 W'(a 2 ) 



we have the same answer as BIZ 



12 



log(2 - a 2 



(4.29) 



(4.30) 



It is important that our result is quite general and we can easily understand the 
generating function for each model using ( ^4.27| ). 

If we adopt W{r) discussed in section 2 



W{r) = 1 + E (j 



r 



T; 



(4.31) 



we can roughly understand the scaling characteristic of the toric generating function. 
Using W(a 2 ) = 1 we have 



k-l 



-W\a 2 ) 



1 - 



k 



E 



a 
~k 



2 \ Mi+2 



1 



X 



Mi + m + 2 - k 
k 



1 - 



a 
' ~k 

k 



rii —k 



(4.32) 



18 



Therefore, 




i=0 




a 

T 



(4.33) 



Indeed e^l) has a logarithmic singularity with respect to 1 — a 2 /k. If one wants to 
see the scaling nature of the generating function more explicitly, a scaling unit like 
the cosmological constant must be introduced. We will argue this point in the next 
section. 



4.2 Double Torus 



Next we derive the generating function of correlation functions for double torus. The 
process of finding the answer in the double toric case gives us a key for generalization 
to higher genus case. 

The starting point for the generating function of double torus is more complex 
than the toric case and is read from ( |1.15| ), 



,(2) 



dx(l — x) 



r± 1 fr 2 \ 2 
r 2 \r J 



192 



22{iy"(0)} 3 - 20W"(0)W"\0) + 3W (4) (0) 

(1) 1 1 r 0) 1 



12 I r 



— {(1 — x) log —1 
720 I V x ) 



(4.34) 



We already know all but r±. Before we express by W{tq), we rewrite last two 
terms of the right hand side of ( [4.34 ). As we see in the toric case we should treat 
r 2 in the form of f 2 = ^2/^0, 



12 l l VnJ 



(i) 



12 L To 
If we use the following relation 



12 r 



1 /r 2 \« 
o~12 



(4.35) 



x=0 



ax ar aro 

and the result for r 2 /r in the previous subsection, we have 



(4.36) 



r 2 \W _ 8{W"(r )} 3 - 7W'(r )W"(r )W"'(r ) + {W'(r )} 2 W^(r ) u ^ 
) " 12{W'(r )} s 



.r 
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Since W(0) = 1, we have 



^ {a-*) 



144 
1 

144 
1 

144 



2{W"(a 2 )} 2 - W'(a 2 )W'"(a 2 )] [W'(a 2 )}- A 

2{W"(0)} 2 - W"'(0) 

8{W"(0)} 3 - 7W"(0)W"'(0) + W w (0)} . 



(4.38) 



On the other hand, we find 



720 



(1 -ar) log — \ 
x ) 



(3) 



720 
1 

720 



„(1)\2 



r (2) 1 , 

Tq X 1 



'o 



r x~ 



W"(r Q ) 



+ 



r 2 {W>(r )}* r {W'(r )} 3 {W(r ) - W(0)} 

3W"(r ) 



(1 W-(n>)) ( r 3{^(r )}3 + rl{W'{r,)Y 
3{W"(r )} 2 - W'(r )W"'(r ) 2 



r {W>(r )}5 



{W(r )-W(0)}\ 



where we used 



{3) _3{W"(r )} 2 -W'(r )W" f (r Q ) 



W'(r Q f 
Taking the limits of r — > a 2 and r — > 0, we find 



(4.39) 



(4.40) 



720 



1 

240 



(1 — x) log 
1 



r \ 



(3) 



iy"(a 2 



2880 
1 

2880 



a 2 {W^'(a 2 )} 3 a 4 {W\a 2 )Y 



9{W"{Q)} 2 -4W"'{Q) 



io{iy"(o)} 3 - 8W"(o)W"(o) + iy (4) (o) 



(4.41) 



Now we rewrite r 4 only by r$ with the help of the function W{tq). The recursion 
relation for r 4 is much longer than that for r^- This relation is found by identifying 
the coefficient of e 4 in both sides of ( |1.8|) , 



= r 4 + ^2(p+l)A ? 

p=i 



2p+l 
P 



^4^0 + r 2 



2p+l 
V 



p-1 
pr 2 r% 
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paths i paths *>J J 

+ r .{( 2p ; 1 )^- i +( 2j '; i )ip(p-i)« 

+(p-i) £ EkwlV + lp-^ £ E'w^J 

paths 1 

+2 4 £ £' 

paths l >i 



(l)^2„p-3 




SiSjT'2 T, 



(1) (1) n-2 



paths *>J 

lij o-a , y- v^I_2 (2) 0-1 

r o + / > / . 2 * 2 
paths * 

paths * ' paths hi 



1 1 V^'^ 2 2/ (2)\2 p-2 , 1 V^'^ 2/ (1)\2 (2) p-3 



2 ^ ^ 4 

paths l J 



2 2 
paths l 'i- fc 



r5 Tq 



4! 



paths i >J\<M 



(4.42) 



In this equation, numbers J p ({a}) given by partitions of two and four appear. For 
ip(2) and / p (l 2 ), we have computed before. Thus we have to determine 5 sums 
characterized by partitions of four: 



/p(2 2 ) 

/ P (i 4 ) 



E E^f. ^- 3 )= E £'^ 3 , 

paths ho 



paths 

I 

paths l .i 



1 1 
_ \ ~> \r^/^„2„2 

O Z^ /I 



2 ^ ^ 4 1 v 



s-s-. i P (i 2 -2) = - E E 



paths » J>fc 



2 



1 

4! 



E E ' S i S 3 S kSl 



(4.43) 



paths ^'.M 

Applying the argument of the previous subsection, we can compute the sum of the 
5 numbers, i.e. , G p (4). Expanding the g-deformed combination 2 P +iC p (x) up to x 4 
order, we have 



paths 



4! 

G P (4) 

J p (4) + J p (l • 3) + J p (2 2 ) + J p (l 2 • 2) + J p (l 4 ) 

{2p+l / p P+l\ 

£ ^ - l)(3r 2 - 7) + ( E + E ) (r 2 - Wr 2 - 3) 



5760 
1 f 1 



2 112 



(P + 1) 2 P 



1 1 

2242 



,r=l r=l, 
'2p+l p p+l\ 

E+E + E (- 2 -i) 2 

r=l r=l r=0 J 



2p+l 
P 
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1 -(p + l) 3 p(5p 2 + 2p - 2) ( 2p + 1 J . (4.44) 



1440 v 7 ' 7 \ P 



Using the generating function r p (x) denned in ( [4.20D , we can also compute I p (4). If 



we expand the g-deformed integers in T p (x) up to x 4 order, we find (see appendix 
B) 

W) = ^E(\ +1 )(2p + l-2n){(2 P +l-2n) 2 -l} 
x{3(2p + 1 -2nf - 7} 
1 r(p+l)p(6p-l) ( 2p+ M . (4.45) 
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How we progress from now on is a problem. The generating functions 2p+\C p {x) 
and T p (x) are insufficient to find all of sums given by partitions of four. We need 
another 3 relations in order to decide all the unknown quantities. It seems that any 
approach using generating functions does not work well. It is possible to find all 
of these sums, however. We can guess beforehand p-dependence of the unknown 
quantities. If one carefully looks at the results for I p {2), / p (l 2 ) and I p (4) and the 
expansion coefficients of the g-deformed combination, G p (2) and G p (4), we observe 
following characteristics of them: 

1. They always contain the combination 2 P +iC p and can be divided by (p + l)p. 

2. G p {2n) is a product of 2 P +iC p and a polynomial of degree 3n. 

3. Each of I p (2), I p (l 2 ) and I p (4) are also a product of the combination 2p+iC p 
and a polynomial of p. The degree of the polynomial for / p (l 2 ) is three, which 
is the same as that for G p (2). It is natural to think that the degree of the 
polynomials are decided by the number of summation indices. Indeed the 
number of indices is both two for G p (2) and / p (l 2 ). On the other hand, the 
degree of the polynomial for I p {2) is less than that for / p (l 2 ) by one, because 
there is only one summation index for I p {2). The degree of the polynomial for 
Jp(4) is also explained by such a consideration. 

4. Since / P (l 2 ) vanishes at p = 1, it must have a factor (p — 1) in the polynomial 
part. In the same way, J p (l-3) and I p {2 2 ) should have (p — 1) in the polynomial 
part. There will be a factor (p — l)(p — 2) in the polynomial part of / p (l 2 • 2) 
and (p - l)(p - 2)(p - 3) for J p (l 4 ). 

Therefore all we have to do is decide two parameters for each of the unknown 
quantities. We assume that 

/ p (l-3) = (p+l)p(p-l)(a 1 p + 6 1 ) 
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J p (2 2 ) = (p+lMp-l)(a 2 p + 6 2 )( 2p + 1 j , 

J p (l 2 -2) = ( p+ l) p ( p -l)( p -2)(a 3 p + b 3 )^ 2p + i y (4.46) 

Unknown parameters a 1; b±, a 2 , 62, o 3 and 63 are fixed by explicit computations of 
I 2 (l ■ 3), J 3 (l • 3) and so on. For example, 

/ 2 (l-3) = ^[2 + 8 + 1 + 1-1-1 + 1 + 1 + 8 + 2] 
11 

y 

/ 3 (l-3) = ^[2{(54 + 24 + 27 + 3 + 8 + 2) + (16 + 16 + 8 + 2 + 8 + 2) 
+2(8 + 2 + 8 + 2 + 1 + 1) + 2(8 + 2) 
+ (8 + 2- 8- 2-1-1) + (1 + 1 + 1 + 1 + 1 + 1) 
+3(1 + 1) + (1 + 1 - 1 - 1 - 1 - 1)} + 4(-l - 1)] 

= f • (4-47, 

It is a laborious work to calculate numbers such as ^(l • 2) by hand. Of course 
these quantities are easily computed if one uses a computer. After a long algebra, 
we find these quantities 

w 2 ) = ^o(p + i)p(p - i)(6p - 1) ( 2p + 1 ) , 

/ P (l 2 • 2) = ^(p + l)p(p - l)(p - 2)(llp - 3) ( 2p + 1 ) . (4.48) 

We checked that these results are valid even for other values of p which we do not 
use to determine the unknown parameters. The rest one has a factor (p + l)p(p — 
— 2)(p — 3) as expected 

/ P (l 4 ) = -^(P + 1)P(P - 1)(P - 2)(p - 3)(5p - 2) ( 2p + 1 ) . (4.49) 

The degree of the polynomial part in I p (l 4 ) is six, which is the same as that of 
Gp(4). These facts support validity of our assumption for the sums. 

Again we can rewrite the recursion equation for r 4 by derivatives of W(r ), 
r 4 W'(r ) + l -W"{r )rl + r 2 ^W"(r )r^ + l^foXr^) 2 } 



+ r 



^W'"(r )r 2 r^ + lw^(r )r 2 (r^) 2 + \w"' {r^ + \w"{r,)rf ] 
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r 1 

1 


1 

T 




.60 r 




1 


i 

+ 




.40 r 




■ 11 


+ 




[360 




■ 1 


+ 




[288 



11 

— 1 

240 



r 2 W^(r ) + -r W^(r ) 

r 2 W^(r ) + ±-r W (5 \r ) 
oU 



rj 4) + 



^r 2 ^( 4 )(r ) + ^r ^"(r ) 



,(2)n2 



0. 



(4.50) 



We rearrange this equation using the relation found in the previous subsection 



lW"(r )r^ + ^-W'"(r )(r 



(1)N2 



6 



12 



--W'Co) 
r 



(4.51) 



and identities 



,(2) 



.(l) r 2 



r o H r — 



r 

.(2) r 2 



r 2 \ 



(i> 



r - 







+ 2r, 



r 



(1) 



+ r 



r J 



(2) 



(4.52) 



and find 



r 



W'(r ) 

Ks) V ^-i(^)W(n,)(r, 

-i (^) W ff (r )r, 



,(2) 




^0 

1 

3 



(1)N2 
) 



r 2 \ 



(i) 



(2)\2 



^ (4) (ro)(r«) 2 rf-liy( 5 )(r )(r«) 



-?'o 



+if^W"(ro) r S 2, + ^('^W (1) M(r< 1) ) 2 



6 Vr 



12 Vr 



+ ^'"(ro)ri 4) + ^ (4) (ro)rfV^ + ^ (4) (r )(rf) 2 



11 



(4.53) 
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At this stage, we have terms proportional to ro and terms which do not depend on 
r explicitly. The relation x = W(r ) and Q4.25| ) found in the toric case allow us to 
write 



f ) = - [l5{iy"(r )} 3 - 10iy'(ro)iy"(ro)iy w (ro) 
+{^'(ro)} 2 ^( 4 )(r )l [^'(ro)]" 7 , 



(4.54) 
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and 



7jA 
To) 



(2) 



^ [48{W"(r )} 4 - 59W'(r ){W"(r )} 2 W"'(r ) 

+7{W'(r )} 2 {W"\r )} 2 + ll{W'(r )} 2 W"(r )W^(r ) 
-{W'(r )} 3 W^(r )} [W'(r )}- S . (4.55) 



We get to the following compact result similar to that for r 2 /r in the end, 

2 



r 4 1 / r 2 x 2 
2 \r , 



(2), 



(4.56) 



V^'(r )rfr 

where functions Q\ (r ) and £1$ (r ) are described only by derivatives of W(r ), 



1 



^ 2) (ro) 
Finally we have 

/ dx(l — x) 
Jo 



t , |Q [28{^"(r )} 3 - 29W'(r )W"(r )W"'(r ) 
+5{W'(r )} 2 W^(r )] [W\r )]' 5 , 

r 31{^"(r )} 2 - lQW'(r )W'"(r )} [W(r )\- 4 . (4.57) 



2880 



r 4 


\ 


©1 


r 







/ dr(l-W(r))- 
Jo i 



(2), 



dr W'(r) dr 

V) + ^ 2) 

a 2 [28{W"V)} 3 - 29W'{a 2 )W"{a 2 )W"'{a 2 ) + 5{W'(a 2 )} 2 W^(a 2 ) 



-a 2 Qf\a 2 ) + Q ( 2 \a 2 ) + fi[ 2) (0) - ^ 2) '(0) - ^ 2) (0) 



1440{W"(a 2 )} 5 
31{W"{a 2 )} 2 - 16W'{a 2 )W'"{a 2 ) 

2880{W(a 2 )} 4 
90{W"(0)} 3 - 84W"(0)W"'(0) + 13g(Q) 
1440 

sijiy^o)} 2 - i6iy w (o) 

2880 



(4.58) 



As in the case of torus, terms that include W"(0), W"'(0) and W^ 4 )(0) miraculously 
cancel out each other, 



l — [90{W^"(0)} 3 - 84W"(0)W"(0) + 13W^ (4) (0)] 
^[31{^"(0)} 2 - 1QW"'(0) 



1440 



2880 



[22{W^"(0)} 3 - 20W"(0)W"'(0) + 3W {4) (0)] 
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1 



144 
1 



+ 



144 
1 



2{w"(o)} 2 -w"'(o) 



i{W"(0)} 3 - 7W"{0)W"'{0) + W^{0) 



2880 
1 

2880 

= 0. 



9{W"(0)} 2 -4W'"(0) 



10{W(0)} 3 - 8W"(0)W"'(0) + W {4 \0) 



Getting all contributions together, we have the generating function for double torus, 



,(2) 



28{W"(a 2 )} 3 - 29W\a 2 )W'\a 2 )W"\a 2 ) + 5{W'(a 2 )} 2 W^(a 2 ) 



1440{W^'(a 2 )} 5 
9{W"(a 2 )} 2 - 4W'(a 2 )W"'(a 2 ) 

2880{W"(a 2 )} 4 
1 1 W"(a 2 ) 1 1 



240 240 a 2 {W'(a 2 )} 3 240 a 4 {W'(a 2 )} 2 

Again we check our result by the example of || 

W(r) =r + \2X l r 2 . 



(4.59) 



Since 
we find 



a 2 W'(a 2 ) 



3 (2) 



2 - a 2 , cfW'ia 2 ) = 2(1 - a 2 ) 



1 28-8(l-a 



2\3 



9 4(1 -a 2 



1440 (2 - a 2 ) 5 2880 (2 - a 2 ) 4 
1 1 2(1 -a 2 ) 1 



240 240 (2 - a 2 ) 3 240(2 - a 2 ) 2 
1 (l-a 2 ) 3 (82 + 21a 2 -3a 4 ) 



720 



(2 - a 2 ) 5 



(4.60) 



This completely agrees with the result of BIZ and guarantees validity of our method 
for computing summations over the paths. If you remember how laborious the 
calculation of BIZ is, you will be satisfied with the convenience of our result (|4.59|). 



4.3 h > 3 Cases 

We succeeded in solving e^ h ' for given W(r ) up to h = 2. It is possible to apply 
our method for computing generation functions to higher genus case. 

The most difficult point is how to compute the sums over paths I p (l ai ■ ■ ■ (2n) a2n ) 
given by a partition of an even integer 2n. These summations are necessary for 
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solving the recursion relation for r 2s defined in ( 1.10| ). Once we overcome this point 



it is straightforward to solve the generating functions for higher genus. Of course, 
it is a hard work from the practical point of view. 

As discussed before, G p (2n) and I p (2n) are both products of the combination 
2p+\C p and a polynomial of p. Our success in double toric case encourages us to 
assume that general I p {l ai ■ ■ ■ (2n) a2n ) has the same structure as G p (2n) or I p (2n), 



■ ■ • (2n) a -) = A(p) , (4.61) 




where A(p) is a polynomial of p. This assumption is reliable. Let us consider the 
degree of the polynomial A{p). As easily confirmed by expanding 2 P +iC p (x) with 
respective to x 2 , the degree of the polynomial part of G p (2n) is 3n. Since the number 
of indices of the summation in G p (2n) is 2n and that in general I p (l ai ■ ■ ■ (2n) a2n ) is 
m defined by ( |4.10| ), the degree of A(p) is 3n — 2n + m = n + m. Furthermore we 
know that / P (l Ql • • • (2n) a2n ) vanishes for p = 1, 2, . . . , m — 1 and A(p) always has a 
factor (p+l)p. Indeed the difference of G p {2n) and I p {2n) has a factor (p— 1), since 
3 Ci(x) = T 1 (x) = N 3 (x). Thus we can assume the detailed form of the polynomial 
part 

I p (l ai ■ ■ ■ (2n) a2n ) 
= A(p)(p+l)p(p-l)---(p-m + l)(^ 2p + i y (4.62) 

where A(p) is a polynomial of degree n — 1. Since this polynomial is determined 
by fixing n coefficients, we need n individual values for I p (l ai ■■■(2n) a2n ), p = 
m, m + 1, . . . , m + n — 1 to determine each of the sum completely. We can express 
r 2h or r 2 h with respect to derivatives of W(r ) using the recursive equation for r 2 h 
and the knowledge for f 2s , (0 < s < h). This means that we can solve for 
any h in principle. In general, the larger the genus is, the harder the computation 
of the sums will become. We have full information for the generating functions in 
exchange for complexity of the calculation. If one wants to compute the generating 
functions for h > 3, one must rely on a computer program to understand the sums 
over the paths. 

Now we guess what happens when one calculates general e^. First, from our 
experience for toric and double toric cases, r2h, (h > 1) will be written in a form 
like 



In (|4.63|) an unknown function i?2/i(^o) will have following structure, 

h-l 



^(ro) = E(- 1 ) S ^f ) (ro), (4.64) 



s=0 
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where f^(ro) are written by only derivatives of W{tq). This is based on the fact 
there is h unfixed coefficients in the polynomial part of I p given by a partition of 2h 
as we showed in (|4.62j ). We have already used these notations in the double toric 



case. For the toric case, 



As we will see in section 5, only O^-i contributes in the double scaling limit. This 
means that only the first term of the right hand side of (|1.15|) contributes in this limit. 
Probably there will be a rule in the form of fi^ (r ) . Our method cannot answer this 
question satisfactorily. Our conjecture for f 2 h is important, for it guarantees that we 



can perform integration in (|1.15 ) explicitly. The case of sphere is the simplest but 
is exceptional, for fo cannot be written in the form of Q4.63Q and there are regular 
terms in e^ -*. 

One of supports of validity of our method is quite a surprising cancellation in 
terms which include derivatives of W(tq) at ro = 0. It is natural to think that these 
cancellation will also happen in the case of h > 3. This conjecture is expected from 
our knowledge that the regular part is present only in genus 0. This cancellation is 
concerned with the fact that there is a function E^ h >{r) such that 



e 



W-EW(a 2 )-E( h \0) 



up to a constant. We can regard — E^ h '(0) as the regular part and E^ h \a 2 ) as the 
singular part. We find such functions for h — 0, 1, 2: 

£?<°>(r) = J^l-W(r)Y, 

nf11/ s 1 , rW'ir) 
E {1 \r) = log 



E^\i 



12 & W(r) ' 

28{W"(r)} 3 - 29W'(r)W"(r)W"'(r) + 5{W'(r)} 2 W^(r] 



lU0{W'(r)Y 



9{W"(r)} 2 - 4W'(r)W'"(r) 
{W'{r)Y 
1 1 W"(r) 



24(W(r) 240 r{W'(r)} 3 240 r 2 {W"(r)} 2 
'bulk' 



(4.65) 



One can easily check that -E^{0) = A,, £«(0) = E^{0) = 0. 



5 Correlation Functions and Leading Singular 
Terms 

In this section we make a study of correlation functions of scaling operators in 1- 
matrix models. It is well understood that these correlation functions obey a scaling 
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rule in the continuum theory ||. Our results for the generating functions of the 
correlation functions are exact. As we do not take any continuum limit, we can 
study the finite size effects of the correlation functions as well as the universality of 
the models. 

We must prepare a scaling unit to understand scaling properties of the correlation 
functions. Here we set up the model by the function W(r) in (|2.9|) . Then the system 
is fine tuned in the k-th criticality and has (m + 1) source terms. We introduce the 
scaling unit as the 0-th scaling operator as announced in section 2. It is convenient 
to choose a puncture operator 

P °° m= 12^ tr$4 (5 ' 1} 

for this purpose. In this section, we use a letter t as the scaling unit instead of r . 
This is nothing but the rescaled cosmo logical constant. The scaling unit t is slightly 
different from the rescaled cosmological constant of t. There is a relation between 
these two scaling units 

t = j\r~ 2fc/(2fc+1) i . 

In the double scaling limit, we take a limit iV —>■ oo keeping t finite. This procedure 
indeed extracts only the universal part. Therefore the most singular contributions 
with respect to t are universal and survive in the double scaling limit. The remainder 
are regarded as finite size corrections. 

We may choose a coupling constant of other scaling operator as a scaling unit. 
However marginal scaling operators whose gravitational dimension are one are not 
suitable for the scaling unit. The reason why marginal scaling operators are not 
accepted is that coupling constants of them cannot scale. We have already met such 
a situation. The fine tuned part of the action S CT is effectively a marginal scaling 
operator and its coupling constant is r_i = —1, i.e. , a number. Such a change of 
the scaling unit brings a seemingly new critical behavior. Effective values of the 
critical exponent and the gravitational dimensions of scaling operators change into 
new values. In the end of this section we argue this point. 

It is convenient to introduce 

/(t,{r}) = l-p W/(a 2 ) = l, (5.2) 

as a function of the cosmological constant t and m coupling constants {r}. The 
condition W(a 2 ) = 1 is translated as a condition for /, 

m hh+2 I w , n \ 

o = E E „ (-W^ 

i=-l p=0 V V J 

= -f + (i - fft + E E „ {-wr**n ■ (s-s) 

i=l p =o \ P J 
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This equation gives a unique solution for f(t,{r}), if one solves ( |5.3|) iteratively. 
Using the technique of the appendix A of 0, the solution of ( |5.3|) is 



f(t,{r}) 



,i/k 



i 



&(i - t l/k rt ni/k - 1 ((73^ - l ) ) ^ 



+t 1/k . (5.4) 

We can understand scaling properties of correlation functions more clearly through 
/. As the first step we rewrite the generating functions for genus 0, 1 and 2 with 
respective to / and its t-derivatives. 

We begin with the case of sphere. Since the bulk part does not contribute to 
scaling properties of correlation functions, only the singular part 

(o) __1 f v tt V ( Mi + Mj + 3\ ( )p +n . +1+p 

(5.5) 

is necessary for us, We can simplify this if we take the second derivative of 
with respect to t. Since t is the cosmological constant, this second derivative stands 
for specific heat. To see this we differentiate once, 

0(0) _1 f Vrr V7_1V/' M * + M i + 3 \ f n l+ n 3 +p f (l) 

dt e smg ~ 2 .2. Z^^Z-A L ) I p J J I 

i=— Xj=— 1 p=0 \ 1 / 

Here we use identities 

^ Mi + Mj + 4 ^ = ^ Mi + Mj + 3 ^ + ^ Mi + Mj + 3 \ ^ ?) 

and (lU(a 2 ) - l) 2 = to derive 



m m 



= ]T £ £(-l) p ( Mi + 5' + 4 ) / ni+nj+p 

i=— 1 j=— 1 p=0 
mm oo 



P 



i=—lj=—l p=0 

Since / is not one, we have 



V v r _ E( _ ir /M i + M, + 3\ /ni+nj . +p(1 _ /) (5 _ 



^ (o) _ _ ^ (-1) P / M + 3 i fni+1+p 



-e 



- E E -^tt^- ' / ni+1+p ■ ( 5 - 9 ) 



at sm § ^ ^ rii + l+p \ V 
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In the same way we differentiate again, 

92 r (°) - _ v r- v 3 (-iv f Mi + 3 ^ p+pfW 
& 2 sing ~ a P =o V p J 



1+p 



P =o 1 +P \P 



m Mi +2 / 

E ^ E (-i) p 

i=-l p=0 V 



f m+P f {1 \l-f) 



P 



~ \[1 ~ (1 - /) 4 ] ■ (5.10) 

Finally we get to 

Ikiig = ~[i - a - /) 4 i = -/ (1 - §/ + / 2 - j/ s ) ■ 

Therefore the function / is almost the specific heat. The deviation of the right hand 
side of ( |5.11| ) from — / represents a finite size effect of the model, since these terms 
disappear in the double scaling limit. Thus is given by 



(0) 

where F(t, {r}) satisfies 



e sing = F (t, M) , (5.12) 



^nt,{r})=-\[i-{i-m- 

Next is the case of torus. In this case, we do not have to differentiate e^. Since 

W(a 2 ) = 1 + E E ( p J ("1) P (l -j) n, (5.13) 
differentiating both sides of ( |5.3| ) once leads to 

= E^ E „ )(-mni+p)f ni+p - 1 f {1) 

i=-l p=0 \ " / 

+ 5(p) ( - 1)p/p 

= -W(a 2 )/ (1) + (1"/) 2 , (5-14) 



where J*™) = d n f/dt n . We have 

/(i) 



kW\a 2 ) = y J> (5.15) 
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and 

k 



a 2 W'(a 2 ) = ^.kW'(a 2 ) 



3 



(1-/) 
/(I) • 

Therefore is written as 

e« = -L\og{a 2 W'(a 2 )} 

= ^log/ (1) -ilog(l-/). 
In the end we rewrite by /. In the same way as we derive kW'(a 2 ) we 

.2„„,^ (l"/) 2 / (2) , 4(1-/) 



k 2 W"(a 2 ) = "—j^r- + 



(/ (l))3 » 

/,- :i ll""(r/ 2 ) ! 



_ (l-/) 2 [3(/( 2 )) 2 -/«/( 3 )] 



(/ (1) ) 5 
| 6(1 ~/)/ (2) | 6 



(/ (l))3 > 
fc 4 w (4) (a2) = (l-/) 2 [15(/ (2) ) 3 -10/ (1 ¥ 2) / (3) + (/ (1) ) 2 / (4) ] 

, 8(l-/)[3(/( 2 >) 2 -/ (1) / (3) ] 12/( 2 ) 

(/ (l))5 + (/ (l))3 • 

If we note = k n W^(a 2 ), 

(2) (1 - /)[28(iy( 2 )) 3 - 29W^W^W^ + 5(iyW) 2 iyW] 

1440(iy( 1 )) 5 
9(VF(2))2 _ 4^(i)|y(3) 



1880(W^)) 4 
1 I^( 2 ) 



240 240(1 -/)(iyW) 3 240(1 - f) 2 (WW) 2 
l 6 (/(2))3 _ 2i/(i)/(2)/(3) + 5(/( 1 )) 2 /W 

1440(/W) 4 (1 - 7p 
2i 5 (/( 2 )) 2 - 156/W/( 3 ) 
2880(/( 1 )) 2 (l -/) 4 
47/( 2 ) 593(/«) 2 1 



120(1 -/) 5 720(1 -/) 6 240' 
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This result is a little complex. We will be able to rewrite the generation functions 
for h > 3 in this manner and the results will be more complicated than the double 
toric case. 



Eqs.( |5.l2p , ( |5.17| ), ( p.21| ) and the formal solution for / (|5.4| ) give us exact answers 



for the correlation functions of scaling operators, 



[Pm/k Pn p /k) - dn q e M=Q . [b.ZZ) 



'p 

Here we do not go into details, for we cannot find any useful things in finite size 
corrections. We will go deep into the leading contributions soon. If one asks for 
exactness, we must study influence of the deformation parameter b discussed in 
section 2. If we define 

Eqs. (|5.12|) , ( |5.17p and ( |5.21|) hold good as they are, while the equation for / is a 



little modified to 

o = -f + ki - f)f k + (i - f?t + EE „ {-iyr +p n . (5.23) 

i=l p=0 v v J 

We can extract leading singular contributions from the exact generating func- 
tions. These contributions survives in the double scaling limit. From now on we 
add a subscript to indicate the most singular contribution and take m = oo. The 
equation for f (t, {r}) can be easily read from ( |5.3| ) by discarding p ^ terms, 

oo 

E fo'n = , (5.24) 

i=-l 

or 

oo 

i=i 

This equation is the same as the one for specific heat given by Gross and Migdal 
(GM) |]J when they studied matrix models on the sphere. Indeed we can easily 
understand that fo is the specific heat, if we pick up the leading contribution from 

dm), 

losing = -fo ■ (5-26) 
Therefore the generating function for sphere is 

<Lg = -W,{t}) , (5-27) 

where F is related to f as (d 2 /dt 2 )F = f . Since our result for eg g- n g is exactly 
the same as that of GM, the correlation functions of scaling operators also coincide 
with their results. We do not repeat any longer. 
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It is important to pay attention to absence of the labels M in ( |5.25| ). This means 
that all representations p^y fc [$] with different M but the same n are degenerate in 
the double scaling limit. This shows that the indices M are non- universal while n 
is universal. Therefore we omit this non- universal index M here and use p n /k [$] for 
the scaling operator with dimension n/k in this approximation. The M-dependence 
is associated with the finite size effect of matrix models. If one wants to study 
finite size effects, one must learn much about the M dependence of the generating 
functions. 

Let us consider an influence of the deformation parameter b which connects a 
criticality with a neighboring one. Formally only the equation for / is modified 
when b ^ 0. According to this modification, the equation for / is modified to a 
form as 

oo 

t=(l-&)/o -E^/o\ (5-28) 



where b ^ 1. This shows that the modified equation (|5.28|) is essentially equivalent 
to the original one ( |5.25|) by simultaneous rescaling of t and {r}. Therefore this 
deformation parameter is non- universal if we do not fine tune it. 

The greatest advantage of our method is that we can compute the leading singular 
terms of the correlation functions in higher genus case. Such correlation functions 
have been formally studied in the double scaling limit 0. It was, however, a hard job 
to compute them explicitly. Our method offers an easier procedure for computing 
the correlation functions for higher genus than the method using the double scaling 
limit. 



From ( |5.17|) and ( |5.21| ) we can write down the leading part of the generating 
functions immediately, 

^log/o (1) , 
16(/ «)» - 21/W/f + 5( /o « )Vo (4) 



e (1) 



(5.29) 



CO 



If one traces terms in the leading contributions to their origin, they come from Q h 
as mentioned in the previous section. Anyway one can pick up the leading terms 
roughly taking / ~ t 1 ' h . As explained in the appendix A of H, we have 



fo(t,{r}) 




P =i 



Thus the derivatives of f are (n > 1 

oo i / a \ P+n-l 



f 
JO 



(») 



k^p\ \dt) 




rii/k 



\i=l 



t 



1/k-l 



m/k \ y.l/fc-1 



t 



l/k 



(5.30) 



(5.31) 
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The leading contribution for the correlation functions of scaling operators are defined 
as 

Ah) _ O O Jh) 



(Pm/k ■ ■ ■ Pn p /k)o 



_d_ 

d7i 



<9t„ 



{t}=0 



(5.32) 



Using ( |5.29| ), ( |5.30| ) and ( |5.31| ), we find the free energies, one-point functions and 
two-point functions for torus and double torus: 



e (1) 



M=o 



1 i , ^-1, 
log k — log t . 

12 s 12k s 



(5.33) 



(Pm/k) 



(1) 



(Pm/kPn 2 /k)o 



(1) 



,(2) 



rii + 1 — k 
12k 



+ni/k- 



12k 2 



[(m + n 2 + 1 - fc)(ni + n 2 + 1 - 2k) 



-(m + 1 - k)(n 2 + l - k)}t {ni+n2)/k - 2 

_ (fc-l)(2fc + 3) 2 _ 1/fc 
{r}=0 720fc 



(5.34) 
(5.35) 



(Pni/fc)o 



(2) 



r (m + 1 - fc)[5(ni + 1 - 2Jfe)(m + 1 - 3Jfe)(ni + 1 - 4Jfe) 

1 — 2A:)(1 — fe)(9 — 2As) 
(5.36) 



1440& 3 

-2l(m + l - 2k)(m + l - 3k)(l -k) + 3(m 
-(l - fc)(n + liife - 2fc 2 )]t ni/fc - 3 - 1/fc , 



(Pm/kPn 2 /k)o 

1 



(2) 



_^(ni+na)/fe-4-l/fe x 



2k)(n 1 +n 2 + 1 - 3k)(m + n 2 + 1 - 4A;) 



2k)(n 1 + n 2 



3fc) 



fc)(9-2Jfe) 



1440& 3 

[5(ni + n 2 + 1 - k)(ni + n 2 + 1 
x (m + n 2 + 1 — 5/c) 
-21 (m + n 2 + 1 - k)(m + n 2 + 1 
x(ni + n 2 + 1 - 4fe)(l - fe) 

+3(ni + n 2 + 1 - fe)(m + n 2 + 1 - 2k){n x + n 2 + 1 - 3k) (1 
-(m + n 2 + 1 - fc)(m + n 2 + 1 - 2k)(l - fc)(ll + llife - 2k 2 ) 
-10(ni + 1 - fc)(ni + 1 - 2fc)(m + 1 - 3Jfe)(ni + 1 - 4A;)(n 2 + 1 - k) 
-21 (m + 1 - fc)(ni + 1 - 2Jfe)(m + 1 
-21(7i! + 1 - fc)(ni + 1 - 2A;)(n 2 + 1 
-10(m + 1 - fc)(n 2 + 1 - k){n 2 + 1 - 
+63(n! + 1 - jfe)(7ii + 1 - 2k)(m + 1 
+96(n! + 1 - jfe)(7ii + 1 - 2k)(n 2 + 1 
+63(m + 1 - £;)(n 2 + 1 - k)(n 2 + 1 - 



- 3k){n 2 + 1 - k){n 2 + 1 
-k){n 2 + l-2k){n 2 + l 
2k)(n 2 + l-3k)(n 2 + l 

- 3k){n 2 + 1 - fc)(l - k) 

- k)(n 2 + 1 - 2Jfe)(l - k) 
2k){n 2 + 1 - 3£;)(1 - k) 



2k) 
3k) 
Ak) 
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-3(m + 1 - k)(m + 1 - 2k){n 2 + 1 - fc)(l - fc)(43 - 22k) 
-3(m + 1 - k)(n 2 + 1 - k)(n 2 + 1 - 2ife)(l - fc)(43 - 22fc) 
+2(m + 1 - k)(n 2 + 1 - fc)(l - fc)(49 - 17 k - 2A; 2 )] . (5.37) 

Especially our results for free energies ( |5.33D and ( |5.35|) are consistent to the asymp- 
totic genus expansion of the whole specific heat given in the appendix C of |3J . If we 
sum up these free energies like (Of), the sum obeys the same nonlinear differential 



equations given in |[|. Moreover we can easily read the leading contribution e ^ 
from the exact generating function e". Once e^ h ' is computed by means of our 
method, we can understand all about correlation functions of scaling operators in 
the double scaling limit. It is important that we do not take any continuum limit 
here. We can reproduce the asymptotic topological expansion of the correlation 
functions in the continuum limit when N = 1. This expansion converges even when 
the cosmological constant t is small, if we take the matrix size N sufficiently large. 

We close this section by arguing other choices for the scaling unit. Let us choose 
a scaling operator p*/ fc , (n ^ k) in place of the puncture operator p{j. Again we 
assign a letter t for the scaling unit instead of r . Then the equation for / is 

OO Nh+2 I yr 9 \ 

o = -f + (i - f) M+2 pt + E E {-ifr^n . (5.38) 

i=i p =o \ P J 
In the same way as we find /-dependence of the specific heat, we have 

§4tg = -/''// 2 "(l-/T M+3 . (5.39) 



Here we discuss only the spherical case for simplicity. When n = M = 0, Q5.39D 
returns to (|5 . 1 1| ) 

To study scaling properties of this system, we pick out leading contributions. 
The equation ( |5.38| ) is simplified to 

oo 

= -/o fc + /o* + E/o n ^. 

4 = 

The solution for fo is 

+t 1/(fc - n) . (5.40) 

This small t expansion is valid only for k > n. When n > k, we must use t~ l as 
a scaling unit. Eq. fl5.40p indicates that the scaling operator p^/k acquires a new 
gravitational dimension 

- _ rij-n _ di - d 
I k — n I 1 1 — do I ' 
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where di = Ui/k, do = n/k. It is important that the dimension of the 0-th scal- 
ing operator and that of the marginal scaling operator are always zero and one 
respectively. The new dimension is allowed to have a negative value. We can find 



another critical exponent like the string susceptibility from the specific heat ( |5.39| ). 
Extracting the most singular contribution from ( |5.39| ), 

x2 f2n+l 

^L e (0) . = (5 41) 

dt 2 °> sm S 2n + l' y ' 

We can define a new exponent 7 in the same way as the usual one, 



&L (o) 

Q t 2 e o,sing 



{r}=0 



oc I t " (f ° r H < k) ' (5 42) 

K 1 V (for n > k) . l °' 4Zj 



From (|5.40|) and fl5.42|) we have 



2n + l 7-2rf 

T = ~Ti 1 = h T\ ■ 5 - 43 

\k — n\ |1 — a | 

Thus the new exponent is always negative as the n = case. The scaling behavior 
of the correlation functions are characterized by the new string susceptibility 7 and 
the new scaling dimensions di of sources. It is natural that there is a universality 
similar to the original one even when we change the scaling unit. 



6 Discussion 



We have solved general 1-matrix models keeping the matrix size iV finite. It is 
wonderful that these random lattice models are exactly solvable. Our explicit an- 
swers for the generating functions offers many information for the discreteness of 
the models as well as the universality. Especially, our method is fit for studying the 
correlation functions of scaling operators. The structure of the bulk terms found 
in section 2 is the most important result which comes from the discreteness of the 
models. In our approach, to take the double scaling limit is equivalent to selecting 
leading singular terms as shown in section 5. Then is it necessary for us to take 
the double scaling limit in order to interpret 1-matrix models as two dimensional 
quantum gravity? In the double scaling limit, all finite size corrections are ignored 
and the general covariance is exactly restored. However, if we require the restora- 
tion of the general covariance only approximately and allow extremely small finite 
size effects, the topological expansion for small cosmological constant is possible 
when the matrix size N is sufficiently large. In this point of view, this approximate 
universality is a consequence of the randomness of the theory. Since the large- N 
expansion is a topological expansion, we can expect the spherical topology is domi- 
nant. This is desirable as a theory of gravity except for approximate restoration of 
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general covariance. Nevertheless we need the results in the double scaling limit, for 
our method is not fit for deriving the hierarchies of the nonlinear differential equa- 
tions for the specific heat function. Moreover an interpretation of matrix models 
as non-critical string theories attracts us. If we can understand the specific heat 
functions of 1-matrix models as systematic as those in the double scaling limit, we 
may not have to take this special continuum limit. 

In our method, there is a mathematically inexact point. Of course, validity of 
our method is guaranteed by many supporting evidences. The point is the assump- 
tion for the p-dependence of sums over the paths (|4.62|) . Though we convince this 
assumption is quite right, a mathematically exact proof is necessary. We must study 
the mathematical characteristics of the sums more deeply. 

Since we solve the recursion relation for r e (|1.8| ), we can find e^ h ' only one by 
one. It is important to integrate our exact results for and find equations like 
nonlinear differential equations for the specific heat function in the double scaling 
limit. It might be impossible, because hierarchies of such differential equations are 
closely related to universality of matrix models. The existence of such equations 
will make our understanding for matrix models rich all the same. 

Though we can solve 1-matrix models exactly, it is nontrivial whether we can 
extend our approach to 2-matrix models. It is very difficult to introduce discrete 
version of the scaling operators systematically even in the Ising case. We have no 
clue to solve this problem. 
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Appendix A g-deformed Combination 



In this appendix we prove Q4.14 ). For the reader's information, we give the first few 
generating functions: 



3 Ci(x) 
7 C 3 (x) 



+ 1 + e x . 



3x , n -2x + 2 e -x + 2 + 2 e x + e 2x + e 3x 



+ e 



e~ 6x + e~ 5x + 2 e" 4x + 3 e~ 3x + 4 e~ 2x + 4 e~ x + 5 



+46* + 4e 2x + 3e* x + 2e* x + e bx + e 



3x 



,6x 
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First of all, we rewrite the definition of the g-deformed combination ( 4.13Q introduc- 
ing y = e^ -1-1 ^, a = s + p + 1, 

p 2p+l 

II = II (l+e"*) 

S=— P (7=1 

2p+l 

= £ 2 p+iC r (*K. (A.l) 

r=0 

One of the new coefficients is also described as a sum over the paths, 

2 P+ A(x)= ]T f[e tr * x = 2p+1 C p (x)#<P +1 >, (A.2) 
paths i=1 

where {at |<7; = 1, 2, . . . , 2p + 1} is another set of p integers and Oi is equal to the 
number of steps until z-th descending down plus one. For the path drawn in Fig.l 
of section 1, they are {a±, cr 2 , cr 3 , cr 4 , <7 5 } = {1, 3, 4, 8, 11}. A relation between Sj and 

Si = Oi -2i, (i = 1,2, ...,p) 
completes the proof of ( |4.14| ), 

e n eSxx = e n e<Tix e- p{p+i)x 

paths i=1 \paths i=1 / 

= 2 P+ iC p {x). (A.3) 



Appendix B Generating Function for I p (2n) 



Here we prove ( [4.20 ) and explain how to compute I p (2) and I p (4). In contrast to the 



proof of ( |4.14| ) in appendix A, the proof of (|4.20|) is diagrammatic. For reference, 
we give the first few of T p (x): 



T^x) = e~ x + l+e x , 

T 2 (x) = e~ 2x + 6 e~ x + 6 + 6 e x + e 2x , 

r 3 (x) = e~ 3x + 8 e~ 2x + 29 e - * + 29 + 29 e x + 8 e 2x + e 3x . 

We explain the first equality in (|4.20|) using Fig. 2. There are (p+ l)p descending 
segments, which are identified by a set of integers indicated in the figure. Let 
us consider a descending segment CD labeled by Since the height of the point 

C is p — i — j, the contributions to T p (x) from the segment CD is the product of 
e (p-*-i)« anc [ tti e number of all paths which pass through the segment CD. On the 
other hand the number of such paths equals to the product of the number of paths 
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from A to C and the number of paths from D to B. These numbers correspond to 
the two binomial coefficients in (|4.20|) . As the generating function T p (x) is the sum 
over all contributions from (p + l)p descending segments, the first equality in Q4.20D 
holds. 




Fig.2 

To prove the second equality in (|4.20|) , we need a device. Comparing the coeffi- 
cients of Q^P~ m ) x in both sides of the last equality of ( |4.2(J| ), the proof resolves itself 
into the proof of 





p — 1 + m — 2i 
m — i 



(B.l) 
(B.2) 

where m = 0, 1, . . . ,p— 1. The case when m = p is special and we must separate as 
( |B.2|) . As we discussed above, the right hand side of ( P-l| ) describes the number of 
the paths which pass through the (m + 1) segments at height p — m marked by a 
circle (see Fig. 3 (a)). Here the paths which pass two marked segments are counted 
twice. In the same way the paths which pass /, (= 1, 2, . . . , m + 1) marked segments 
are counted I times. We turn back each of such paths on the line of height p — m 
(the broken line of Fig. 3 (a)) fixing the part from the point A to the point such that 
the path cuts the broken line first. 
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p-m 




-1 



height 




(a) height 

2p-m - 
2p-2m - 
p-m - 



(b) 




Fig.3: A path which passes through at least one of the marked segments is 
described by a thick line in Fig. 3(a). This path is 'folded' on the broken line 
of height p — m and turns to the path described by a thick line in Fig. 3(b). 
After this folding procedure, all (to + 1) marked segments locate just above 
the broken line. The marked segments which the path passes does not change 
before and after this folding procedure. In this example, p = 4 and m = 2. 

Each of the 'folded' path corresponds to each of the paths which connects A with 
B' as illustrated in Fig.3(b). The marked segments are transformed into the marked 
ones in Fig.3(b). Since the number of paths from A to B' is 2p+iC m , the number of 
the paths in Fig.3 (a) which pass through the marked segment at least once is also 
2p+\C m - This is nothing but the contribution in the left hand side of (|B.1|) from 
n = m. If one notices that there is a overcounting in the right hand side of ( B.l| ), 
it is easy to guess that the number of the paths which pass through the marked 
segment at least I times is 2 P +iC m +i-i- If this is true, we can prove (P-l| ) taking 
n = m + 1 — I, (n = 0, 1, . . . , m). This guess is easily justified by folding each of 
paths in the same way as before moreover 2(m — n) —2(1 — 1) times shifting up the 
height of line for folding by one. Each of the resultant paths corresponds to each of 
the paths which connects diagonal points of a lattice of size n x (2p + 1 — n). We 
can prove (B.2|) in a similar way. In this case, each of paths which passes through 
the marked segments at least I times must be turned back 21 times (/ = 1,2, ... ,p) 
shifting the height of line for folding by one. 

Next we compute I p (2), 

p-i 



Jp(2) 




(2p+ 1 - 2n){(2p 



2n) 



!}■ 



We put P = 2p + 1 and fix polynomials a(P), b(P), c(P) and d(P) defined by 
(P - 2nf - (P - 2n) = a(P){P - n)(P - 1 - n)(P - 2 - n) 
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+b(P)(P - n)(P - l-n) 
+c(P)(P - n) + d(P) . 

If we substitute n = P, we get 

d(P) = -P(P 2 - 1) . 

In similar way, we find 

c(P) = Q(P -l) 2 , b{P) = -12(P-2), a(P)=8, 
If one uses identities like 

/ 2 P + r 



(2p + 1 -n)(2p-n)(2p- 1 - n) = (2p + l)(2p)(2p - 1) 



the quantity of interest is 



i_(2) = — 
P \ ; 24 



P-i 



2p-2 



(2p+l)(2p)(2p-l)£ „ 

ra=0 \ 

-12(2p-l)(2p+l)(2p)X:( 2 ^ 1 

n=0 \ 

+6(2p) 2 (2p+l)X:( 2 r f ) 
-(2p + 2)(2 P + l)(2 P )£^ 2p + 1 ) 

4(2p+l)(2p)(2p-l) |2 2p - 2 + ^ 
-6(2p - l)(2p + l)(2p) x 2 2p ~ x 
+3(2p+l)(2p) 2 |2 2 --( 2 j)} 

-J(2p + 2)(2p + l)(2p) U 2 ^ 1 - 2 f 2p + 1 



1 

24 



2p-2 
p — 1 



^(2p+l)(2p) 



4(2p-l) 



2p-2 
p — 1 



P 



— 6p 



(B-3) 



2p-2 
n 



2p 
P 



+ (2p + 2) 



2p + 1 
p 



(B.4) 



It is interesting that all terms which have a power of two cancel each other. Thus 
we have 



(B.5) 
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The computation of I p (4) is performed in the same way. The cancellation appeared in 
ip(2) case also occurs when one computes L p (4). It seems to a common characteristic 
in computing I p (2s). Of course we can calculate L p (2s) by the parameter fitting 
discussed in §4.3. Probably this parameter fitting will be easier than the method of 
this appendix. 
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